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We examined the diversity and potential drivers of dragonfly distribution in a biodiversity hotspot of 
Southern Africa (Eastern Highlands, Zimbabwe) by surveying 30 sites (13 lentic and 17 lotic habi- 
tats) located within this region. Additionally, we identified the anthropogenic factors that may threaten 
Odonata diversity and abundance. Our results revealed that 27 odonate species are associated with dams 
and ponds, one species is associated with streams and four species are associated with swamp forests. 
Considering odonate diversity between protected and unprotected areas, we found significantly higher 
diversity of odonates in streams in protected areas compared to non-protected areas. Broadly, we found 
anthropogenic activities (e.g. commercial exotic tree plantations, mining activities, vegetation clearing) 
possibly affect Odonata diversity (by decreasing abundance and diversity) compared to those within pro- 
tected and undisturbed habitats. Our results reveal that several human activities like human encroachment 
into riverine habitats potentially threaten the existence of freshwater species in this biodiversity hotspot 
and biodiversity hotspots elsewhere. 
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Introduction 


The Eastern Highlands Biodiversity Hotspot is a hub of endemism in Zimbabwe; known for 
its high number of endemic species of fauna and flora (Branch, Bayliss, & Tolley, 2014; Dalu 
et al., 2016). The highlands are of high economic significance to Zimbabwe due to the number 
of rivers and tributaries that originate from the mountains (Jonnalagadda & Mhere, 2001; Terink, 
Droogers, & en Aa’s, 2014). These rivers in Zimbabwe, congruent to rivers in the rest of the 
world, play a vital role in human livelihoods (sensu Darwall, Smith, Tweddle, & Skelton, 2009; 
Lamsal, Pant, Kumar, & Atreya, 2015) and are undergoing substantial alteration due to climate 
and land use changes from an increasing human population (Zuo et al., 2016; Mukheibir, 2017). 
Mitigating the pressures facing freshwater systems in the world requires scientists and environ- 
mental managers to understand the ecological mechanisms, and ultimately the factors affecting 
biodiversity in ecosystems (Lopes-Lima et al., 2017; Sterling et al., 2017). 
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Odonata species have been identified as one of the priority taxa for environmental health 
and conservation (sensu Clausnitzer, Simaika, Samways, & Daniel, 2017; Hemalatha, 2017). 
Odonates are particularly sensitive to human disturbance and their diversity as larvae and 
adults tends to decrease whenever there are high levels of human disturbances and pollution, 
thereby making them suitable quality indicators for both terrestrial and aquatic ecosystems 
(de Oliveira-Junior et al., 2017; Kemabonta, Adu, & Ohadiwe, 2017). The top-down effects 
(predation capacities) of odonates have also enabled them to be used as biological control 
agents for disease-spreading pests like mosquitoes (Bhusnar & Sathe, 2017; Faithpraise et al., 
2014). 

However, Odonata species are not widely used in Zimbabwe as habitat quality indicators or 
as biocontrol agents because many of the endemic species are yet to be described; with many 
areas of Zimbabwe still poorly explored (Darwall et al., 2009; Dijkstra & Clausnitzer, 2014; 
Dijkstra, Clausnitzer, Méziére, Kipping, & Schiitte, 2011; Samways, 2002; Suhling, Samways, 
Simaika, & Kipping, 2009). There are limited records on the abundances and species distribu- 
tions of most species from Zimbabwe, with most records dating back at least 30* years (Dijkstra 
& Clausnitzer, 2014; Kipping, Dijkstra, Clausnitzer, Suhling, & Schiitte, 2009). The scanty data 
of odonate species descriptions in Zimbabwe may very well be true for many other African 
countries (Bouchelouche, Kherbouche-Abrous, Mebarki, Arab, & Samraoui, 2015; Dumont & 
Verschuren, 2005; Dijkstra et al., 2011) as well as across the globe (e.g. Dijkstra & Clausnitzer, 
2014; Khan, 2015). Consequently, Odonata species have increasingly acquired a high conserva- 
tion status and are commonly included in red-lists of many countries (Dijkstra & Clausnitzer, 
2014; Samways, 2002). 

To date, more than 120 species have been recorded, from more than 170 localities in the 
Eastern Highlands Biodiversity Hotspot of Zimbabwe (Kipping et al., 2009; Pinhey, 1984). 
Furthermore, two endemic species to this region, Africallagma cuneistigma and Pseudagrion 
vumbaense, are hardly studied yet they carry a high conservation importance; because knowl- 
edge of these two species remains rudimentary (Darwall et al., 2009; Dijkstra & Clausnitzer, 
2014). Africallagma cuneistigma was first recorded in 1958 in the Bundi river (Pinhey, 1961), 
Chimanimani mountain and was again recorded in 2013 at the same locality (Dijkstra & Kipping 
2013). Pseudagrion vumbaense was first recorded in 1958 in Vumba and again in 1966 in the 
Chirinda forest in the Eastern Highlands of Zimbabwe. It was more recently recorded in 2013 at 
four localities around the biodiversity hotspot (Dijkstra & Kipping, 2013). 

Most Odonata species recorded in the Eastern Highlands of Zimbabwe were recorded from 
lotic habitats (Clausnitzer et al., 2012; Kipping et al., 2009) with very few studies considering 
species diversity in lacustrine environments. While these aforementioned surveys have con- 
tributed to our knowledge of species distributions in Africa, in most cases the data lack pertinent 
information on species ranges, habitat requirements, and threats to the species (e.g. land use 
practices); aspects that are crucial to the conservation of biodiversity (Allan & Castillo, 2007). 
Furthermore, a few studies have been conducted in Zimbabwe to assess the effect of industrial 
and sewage effluent on Odonata (sensu Dube, Makaka, & Sibanda, 2010). Elsewhere, thermal 
tolerances of Odonata larvae families have also been documented (Lake Kariba, Zimbabwe; 
Makaure et al., 2015). A few other ecological studies on Odonata species have been conducted 
in other African countries (e.g. Damm, Dijkstra & Hadrys, 2010; Hart, Bowker, Tarboton, & 
Downs, 2014). Understanding the factors that influence the distribution and species responses 
to perturbations is necessary as it helps ecologists to prioritize the best management practices 
for sustaining habitat quality, biodiversity, and ecosystem services within biodiversity hotspots 
(sensu Allan & Castillo, 2007). 

The high endemism in the Eastern Highlands Biodiversity Hotspot has not been spared by 
anthropogenic activities associated with the increasing human populations and human encroach- 
ment into protected areas. Zimbabwe’s biodiversity hotspot is under intense pressure from 
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agriculture, illegal logging, and small-scale mining activities, which pose a threat to both ter- 
restrial and freshwater habitats and the diverse species therein (Gandiwa & Gandiwa, 2012). For 
instance, when gold was discovered in the alluvial river sands in the upland grasslands of the 
Chimanimani mountain in the year 2000, over 10,000 small-scale miners from Zimbabwe and 
Mozambique invaded the protected areas and began digging in and around riverbeds. These 
mining activities have led to the destruction of the riparian vegetation, while discharge and 
depth of the rivers have been altered as evidenced by aerial images analysed by researchers 
(e.g. Ndunguru, Dondeyne, & Mulaboa, 2006) 

Presently, the endemic species found in the Eastern Highlands of Zimbabwe are listed as 
threatened (Dijkstra & Clausnitzer, 2014; IUCN, 2019. The endemic Odonata species and other 
freshwater species face potential threats due to the destruction of their habitat by the anthro- 
pogenic induced activities (Clausnitzer et al., 2009; Harisha & Hosetti, 2017). The IUCN (2019) 
asserts the need for more data to be collected on general species distribution, which will urgently 
confirm species distribution range and population status for many species that have been recorded 
in the eastern highlands, including the two endemic species Africallagma cuneistigma and 
Pseudagrion vumbaense. To this end, we addressed the following aims: (i) to describe Odonata 
species found in the Afromontane region of Zimbabwe in terms of total abundance and diversity; 
(ii) to determine habitat associations of the Odonata species in the Eastern Highlands Biodiver- 
sity of Zimbabwe; and (iii) to identify and assess the natural and anthropogenic threats facing 
these and other Odonata species. 


Material and methods 


Study area 


The study was conducted in the Eastern Highlands of Zimbabwe; a biodiversity hotspot along 
the Chimanimani—Inyanga Biodiversity Corridor (Figure 1). The study sites were located in five 
protected areas: Nyanga National Park, Chimanimani National park, Vumba Botanical Garden, 
Chirinda Forest Botanical Reserve, the Stapleford Forest along with adjacent communal areas 
along the corridor. All sampled habitats were characterized and marked during a preliminary 
survey in August 2016. Following this survey, habitats were mapped and grouped into lotic 
(seeps, springs, creeks, streams and rivers) and lentic habitats (dams, pools, marshes, lakes, bogs, 
etc.) according to the physical characteristics of the waterbody, and the availability of water 
resources. Thirty sites were selected, comprising 13 lentic and 17 lotic habitats (Table S1 in 
Supplementary File). Surveys were conducted between October 2016 and June 2017 with each 
site sampled once in each of the three seasons recognized in Zimbabwe; i.e. in October 2016 (hot 
and dry season), January (hot and wet season) and June 2017 (cold and dry season); depending 
on the accessibility of the sites. 


Field collections 


Standardized, timed (30-min) surveys were used to survey adult odonates. All sampling was 
conducted between 1000 hours and 1600 hours, when Odonata activity is high. The surveys 
were conducted by two observers (two experienced Odonata taxonomists) and one recorder, 
walking through all the riparian habitats (intensively searching) and recording individuals flying 
over water and emergent vegetation, or within 2 m of the water’s edge. A 5 x 500 m transect was 
used at each site and the Odonata species encountered were identified and counted. Identification 
of species was predominantly done for each species using close-focusing binoculars and field 
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Figure 1. Map of the study area (Eastern Highlands of Zimbabwe) showing sampling points (1-30) as listed in 
Supplementary Table 1. 


guides by the two observers to avoid misidentifications. Where necessary, one individual of each 
species was caught with sweep nets and examined using a hand-lens to confirm identification and 
subsequently released. Furthermore, species were photographed to provide a permanent record 
of identification and occurrence. Where further identification was required, voucher specimens 
were captured and transported in labelled paper envelopes for further identification using regional 
keys. 


Data analysis 
Diversity and abundance (Objective 1) 


To address objective 1, the Shannon diversity index (H) was used to provide a measure of relative 
diversity. Diversity indices between habitat types and amongst elevation zones were compared 
using Mann-Whitney and Kruskal-Wallis tests respectively. To address the potential signifi- 
cance of protected areas for freshwater species, we compared diversity between protected and 
unprotected areas. 


Associations between Odonata species and biotopes (Objective 2) 


Diversity of Odonata is known to change with elevation (Samways, 2008); therefore, we grouped 
sites into elevation zones according to Timberlake (1994a) based on the elevation of the sites, 
in metres [high (> 1800m), high to medium (1400—1800m), medium (1000-1400 m) and low 
(< 1000 m)]. Subsequently, we used a Kruskal—Wallis test to compare Odonata diversity across 
elevation zones. Furthermore, we grouped sites into biotopes [swamp forests, small mountain 
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streams with calm pools, fluctuating pools and marshes, permanent pools, lakes, vleis or marshes, 
sluggish rivers with quiet pools, artificial ponds and dams, swift rivers and streams with invasive 
alien trees] (Samways, 2008; Table S1). We used a multivariate approach [principal component 
analysis (PCA)] to determine the association between the Odonata species and sampled biotopes. 
To further augment the results from the PCA analysis, we also used an indicator value (IndVal) 
analysis (Dufréne & Legendre, 1997) to assess the indicator species for different habitats. IndVal 
analysis has been advocated as a useful method for assessing habitat associations between species 
and the habitats they rely on (Podani & Csányi, 2010). The significance of IndVal indices was 
assessed using 10,000 Monte Carlo permutations. 


Identifying and assessing the anthropogenic threats (Objective 3) 


The anthropogenic factors were identified during the preliminary surveys (habitat characteriza- 
tion). Transects were therefore also set up at sites that fell within: (i) protected areas reforested 
using exotic trees; (ii) communal areas; (iii) large scale exotic tree plantations; (iv) illegal mining 
areas; and (v) protected areas where vegetation is regularly cleared and modified for tourism. 
Subsequently, to assess the effects of anthropogenic factors on Odonata diversity, we com- 
pared these sites against the diversities recorded from the protected and undisturbed sites using 
Mann-Whitney tests. 


Results 


General diversity and abundance of Odonata (objective 1) 


From the 30 sites sampled, we recorded 81 species (27 Zygoptera and 54 Anisoptera) and 10 
families from 2678 individual adults (Table S2 in Supplementary File). Family Libellulidae 
was the most dominant species (44 species), with the Coenagrionidae family (18 species) being 
the second most dominant. The most ubiquitous species was Trithemis arteriosa (374 individ- 
uals) followed by Pseudagrion kersteni (319 individuals). Only one individual of Zygonoides 
fuelleborni was recorded along the Nyamukwarara River. 

Within protected areas, lotic habitats had higher odonate diversity (H = 2.13) compared to 
lentic habitats (H = 1.76; U=63; Z =-2.02; p <0.05). However, in unprotected areas, lentic 
habitats had higher diversity (H=2.10) compared to lotic habitats (H = 1.60; U = 123; Z=- 
2.94; p < 0.05). Further, lotic habitats in protected areas had higher odonate diversity than outside 
the protected area (U = 85; Z = -2.81; p < 0.05), whereas lentic habitats within protected areas 
had lower odonate diversity compared to those found outside the protected area (U = 75.5; Z = — 
2.53; p < 0.05; Figure 2). 

Considering odonate diversity at different elevations, significant differences (Kruskal—Wallis 
H=8.3, p<0.05 df=3, N=75) in Odonata diversity occurred between medium and low 
elevations zones (Figure 3). 


Habitat associations (objective 2) 


Anisoptera species were ubiquitous (PCA; explains 47.3% of the total variation) with no dis- 
cernible biotope associations (Figure 4a). The majority of Zygoptera species were associated 
with artificial dams and ponds (PCA; explains 66.8% of the total variation; Figure 4b) includ- 
ing Ischnura senegalensis, Azuragrion nigridorsum, Ceriagrion glabrum, Africallagma fractum, 
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Figure 2. Odonata species diversity within and outside protected areas (per each habitat type). Different letters denote 
significant differences (p < 0.05) among habitat comparisons. Error bars are based on standard deviations. 
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Figure 3. Odonata diversity across the elevation zones [high (> 1800m), high to medium (1400-1800m), medium 
(1000-1400 m) and low ( < 1000 m)]. Different letters denote significant differences (p < 0.05) among elevations. Error 
bars are based on standard deviations. 


Phaon iridipennis, Agriocnemis exilis and Pseudagrion glaucum. Platycypha inyangae, Chloro- 
cypha consueta and Pseudagrion salisburyense were associated with small mountain streams 
with quiet pools and streams with invasive alien trees. Pseudagrion hamoni, Pseudagrion hageni, 
Platycypha caligata, Allocnemis marshalli, Lestes virgatus, and Lestes plagiatus were associated 
with sluggish rivers with quiet pools and swift rivers. Elattoneura glauca was associated with 
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Figure 4. Biplots based on a PCA illustrating the associations between the aquatic biotopes (adapted from Samways 
et al., 2005) and (a) Anisoptera (first two axes explain 47.3% of the total variation) and (b) Zygoptera species (first two 
axes explain 66.8% of the total variation) recorded from the Afromontane region of Zimbabwe between August 2016 
and June 2017. 


streams with alien invasive trees. Pseudagrion vumbaense, Pseudagrion kersteni, Pseudagrion 
spernatum and Africallagma glaucum were found to not have any specific biotope requirements 
(Figure 4b). 

In support of PCA analyses, IndVal analysis revealed that Ischnura senegalensis was asso- 
ciated with streams with invasive alien trees, while Proischnura subfurcata, Africallagma 


glaucum, Pseudagrion assegaii and Agriocnemis exilis were associated with swamp forests 
(Table 1). 
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Table 1. IndVal analysis of species that associated with specific habitat types. Only significant 
(p < 0.05) habitat associations are presented. The significance of IndVal indices was assessed using 
10,000 Monte Carlo permutations. 


Species Habitat Indval p Freq 





Lotic 

Ischnura senegalensis Streams with invasive trees 0.42 0.029 4 
Proischnura subfurcata Swamp forests 0.70 0.008 4 
Africallagma glaucum Swamp forests 0.58 0.022 7 
Pseudagrion assegaii Swamp forests 0.50 0.048 1 
Agriocnemis exilis Swamp forests 0.42 0.033 2 
Lentic 

Ischnura senegalensis Permanent pools, lakes, 0.89 0.001 12 


vleis or marshes 


Anthropogenic threats (objective 3) 


Lotic habitats located within the protected and undisturbed sites had a significantly higher diver- 
sity compared to the sites located within protected areas, where vegetation was regularly cleared 
and modified for tourism. The protected and undisturbed sites had significantly higher diversity 
compared to sites located in the communal areas, in plantations with extensive mining activ- 
ities. The endemic species, P. inyangae was recorded at a site that fell outside the protected 
area along the Haroni River. Communal settlements and related activities including fishing using 
mosquito nets and stream bank cultivation were identified at this site. Another endemic species, 
Pseudagrion vumbaense, was also recorded outside the protected area in swampy area near Chik- 
ware dam. However, it is worth noting that this area had been reforested with exotic pine trees 
(K. Mafuwe, personal observations, June 2017). 


Discussion 


Diversity and distribution (Objective 1) 


The Libellulidae family had the greatest number of species, supporting the notion that this 
family is generalist, widespread and can easily migrate and adapt to a wide range of habitats 
(Van Huyssteen & Samways, 2009). Libellulidae as a group are known to be resilient to chang- 
ing environmental conditions (Van Huyssteen & Samways, 2009). Trithemis arteriosa was the 
most abundant species (Table S2) probably because this species can inhabit a range of habitats, 
including both open and shaded forests along lotic systems (Damm & Hadrys, 2012). Regard- 
ing Zygoptera species, Pseudagrion kersteni was the most abundant species, in accordance with 
IUCN reports that suggest that Pseudagrion kersteni is widely distributed across Africa with no 
documented threats (IUCN, 2019). 

Of the four endemic species recorded in the study region (Dijkstra & Clausnitzer, 2014), 
only two were found during this study: Platycypha inyangae and Pseudagrion vumbaense. 
Africallagma cuneistigma and Elattoneura lapidaria were not recorded during sampling, prob- 
ably due to declining populations. However, an increased sampling effort is needed to test this 
notion. Because of the rarity of Elattoneura lapidaria, which is known to occur along open rocky 
streams and patches far from streams (Dijkstra & Clausnitzer, 2014), more targeted sampling is 
needed to find this species. 

We found odonate generalists and heliophilic species to be ubiquitous (e.g. T. arteriosa and 
P. kersteni) along the biodiversity hotspot. The dominance of generalists may indicate land use 
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changes from agriculture that favour generalist species than natural habitats that favour special- 
ized species (Devictor, Julliard, & Jiguet, 2008). Researchers have documented T. arteriosa to be 
most abundant in agricultural habitats (e.g. Seidu, Danquah, Nsor, Amaning Kwarteng, & Lan- 
caster, 2017). We only recorded a single specimen of Z. fuelleborni at the Nyamukwarara river 
(the widest river in our study). Zygoniodes fuelleborni, as a specialist species that inhabits wide 
rivers (Dijkstra & Clausnitzer, 2014), was expected to be present in this Nyamukwarara river. 
However, the chemicals from illegal mining activities that discharge chemicals into the Nya- 
mukwara river (K. Mafuwe, personal observation, June 2017; Chipangura, 2017) may explain 
the low abundance of Z. fuelleborni recorded in this study. 

In protected areas, lotic systems had a significantly higher diversity of odonates compared to 
lentic habitats (Figure 2). Dams and pools in protected areas of Zimbabwe are used mainly for 
tourist water supply and recreational activities (e.g. angling and swimming). Therefore, there 
are more anthropogenic activities around these lentic habitats that possibly alter natural habitats 
and behaviour of odonates. In contrast, the lentic habitats in the unprotected areas had a sig- 
nificantly higher diversity of odonates compared to the lotic systems. This may be attributed to 
the increased supply of artificial water bodies built in several communal areas, which may subse- 
quently increase local dragonfly diversity (Deacon, Samways & Pryke, 2018; Simaika, Samways 
& Frenzel, 2016). 

Pools, marshes and lakes in unprotected areas had a significantly higher species diversity com- 
pared to species found in similar habitats in the protected areas; this may be explained by the 
higher concentration of lentic water bodies, particularly dams in Zimbabwe found in unprotected 
areas compared to protected areas. For instance, Van der Mheen (1994) reported only 87 dams out 
of 9818 to be found in national parks and forests compared to the other land use classifications 
in Zimbabwe. 

Despite having fewer sampling sites of lotic habitats in protected areas compared to unpro- 
tected areas, the species diversity of odonatans found at lotic habitats in protected areas were 
significantly higher compared to similar habitats outside protected areas. This may be due to 
the impacts on the quality of rivers by the increased supply of artificial water bodies such 
as dams as they tend to alter aquatic ecology and river hydrology upstream and downstream, 
thereby reducing the quantity and breeding grounds for a number of aquatic species including 
Odonata, along the lotic habitats (Helland-Hansen, Holtedahl, & Lye, 1995). Additionally, the 
low diversity of odonatans along the rivers in unprotected areas may also be attributed to the 
high levels of stream bank cultivations that is apparent in many areas in the Eastern Highlands, 
such as the Pungwe River basin, as well as the effluent from illegal mining activities that occur 
along rivers such as the Nyamukwarara River and Mutare River (Chipangura, 2017; Nyikadzino, 
2016). 

Odonata diversity at low elevations was greater than at high to medium elevation zones. These 
differences in elevation could be explained by the fact that the forests within the high to medium 
and the mid-elevation zones have fertile soils that are extensively used by subsistence farmers; 
vegetation is constantly being cleared to accommodate these farmers (Timberlake, 2002). Addi- 
tionally, Timberlake (1994a, 1994b) documented that poaching, collection of firewood and the 
gathering of traditional herbal medicines that occurs at mid-elevation forests is major threat to 
biodiversity in these zones. The Chirinda forests were previously much more widespread, hav- 
ing been cleared for agriculture over the last 100 years (Mujuru & Kundhlande, 2007). Further, 
higher species diversity of Odonata observed in the low-elevation zones could be attributed to the 
extreme temperatures and weather conditions experienced in low-elevation zones, that include 
lower rainfalls and higher temperatures often exceeding 35°C in summer (Childes & Mundy, 
1998). These climatic conditions restrict farming in low elevation zones and they consequently 
have experienced fewer human disturbances in the past; this has helped to maintain the thick 
forests within these low-elevation valleys. However, there are still patches of unprotected forest 
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in these areas that could be lost to communal settlements. Human expansion may soon have an 
impact on the biodiversity within the valleys (Gandiwa, Matsvayi, Ngwenya, & Gandiwa, 2011). 


Habitat associations (Objective 2) 


Anisoptera species were widespread across the region and occurred within most biotopes. This 
could be explained by the fact Anisoptera, particularly the Libellulidae family and Aeshnidae, 
contains a great number of species, e.g. Pantala flavescens, that are widespread and general- 
ist species known to occur across a wide spectrum of habitats and microhabitats (Haritonov & 
Popova, 2011; May 2013). 

Our data revealed that most Zygoptera species were associated with artificial dams and ponds. 
Damselflies are smaller species that do not migrate for long distances and therefore prefer 
perennial habitats that have a long-term supply of water and sites for perching (Rouquette & 
Thompson, 2007). Artificial ponds and dams act as reservoirs for rainwater, especially during dry 
periods (de Santos Loureiro, 2012; Kosterin, 2010). As such, they provide microhabitats, perch 
sites and artificial breeding sites that can support many species; this is in agreement with other 
researchers (Oertli et al., 2002). Moreover, artificial ponds and dams were the most abundant 
freshwater habitats identified and surveyed during this study (Table S1). 

Ischnura senegalensis was associated with streams with invasive alien trees. Exotic trees, 
especially wattle, pine and eucalyptus trees, are common in rivers, due the exotic timber plan- 
tations found in the Eastern Highlands of Zimbabwe (Tembani, Madhibha, Marunda, & Gapare, 
2014). These habitats are characterized by extensive human disturbance from logging by tim- 
ber companies and illegal settlers that have invaded parts of the commercial forests (Jimu & 
Nyakudya, 2018; Zisadza-Gandiwa et al., 2014). These conditions tend to make these rivers 
depauperate in dragonfly species, and even threatening some with extinction (Samways, 2008; 
Samways & Grant, 2006). However, I. senegalensis is very widely distributed around the globe, 
and it is known to be highly tolerant to disturbed and polluted conditions and may be present 
in trampled waterholes (Lambertz & Schmied, 2011; Samways, 2008), which could explain its 
association with these habitats. 

Proischnura subfurcata, Africallagma glaucum, Pseudagrion assegaii and Agriocnemis exilis 
were associated with swamp forests in accordance with the habitats described by Dijkstra & 
Clausnitzer (2014), who reported all species to occur in grassy, marshy and swampy habitats 
that are often stagnant and temporary. Proischnura subfurcata was reported to particularly occur 
in marshes found in mountainous areas. Samways (2008) did, however, report two specimens 
of P. subfurcata to occur in grassy and reedy reaches of streams found in non-mountainous and 
low-lying areas of South Africa, though the need for confirmation was noted as it was uncertain 
that these two records were indeed correct identification of the species. 


IUCN evaluations and threats to Odonata 


Three species that are considered vulnerable according to IUCN (2019) were recorded in the 
Eastern Highlands of Zimbabwe during this research (Table S2). We did not find the following 
endemic species: Africallagma cuneistigma and Elattoneura lapidaria. Further sampling (includ- 
ing active searches for naiads) is needed to determine the distribution and range of these species. 
This information will be essential to updating the IUCN Red List on the status of these species; 
this is particularly true considering that the Eastern Highlands will continue to experience huge 
land use changes related to a continuously increasing population. 

Commercial tree plantations, mining activities, vegetation clearing and landscape modifica- 
tion for tourism and communal settlements were identified as potential threats to odonates in the 
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Eastern Highlands Biodiversity Hotspot of Zimbabwe (K. Mafuwe, personal observation, June 
2017). The impacts appear to be greatest in areas whose trees have been harvested and cleared. 
These anthropogenic factors can potentially threaten the existence of freshwater species includ- 
ing the endemic species, P. inyangae, that we recorded at a site that fell outside the protected 
area along the Haroni River. Seidu et al. (2017) documented that disturbance causes reductions 
in stenotopic and perching odonate species assemblages due to the loss of the complex struc- 
ture of the forest canopy. Our findings and those of other researchers (e.g. Biervliet, Wisniewski, 
Daniels, & Vonesh, 2009) suggest that deforestation does not only render habitats unsuitable for 
the colonization and reproduction for forest specialist species, but also threatens forest freshwater 
biota. 


Summary and additional considerations for conservation 


Our results illustrated how the diversity and distribution of Odonata species differs among habi- 
tats with varying land use activities. The findings have direct significance for our understanding 
of the environmental pressures on freshwater diversity, and broader implications for a region 
that is considered the source of endemism in Southern Africa. We recommend the inclusion of 
several diversity studies in different landscape characteristics to unravel some of the other land 
use changes that can affect freshwater animals. 
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